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ABSTRACT 

^ Theories  of  rough  surface  scatter  and  gas  bubble  behavior  are  used 
with  the  Pierson-Moskowitz  wind-wave  spectrum  and  an  empirically-guided 
formulation  of  bubble  concentrations  at  sea  to  calculate  the  true  velocity 
gradient  and  losses  *at  the  surface'^.  These  values  are  then  entered  into 
Bucker's  wave  theory  solution  for  sound  propagation  and  leakage  in  a 
surface  duct.  Curves  of  propagation  loss  are  calculated  for  corparison 
with  ocean  test  data  obtained  with  the  SQS-26  sonar.  The  predictions  are 
shewn  to  be  significantly  better  than  those  based  cm  the  enpirical  equations 
of  project  AMDS. 


I.  INTRODUCTION 


The  purpose  of  this  paper  is  to  show  that  physical  descriptions  of 
bubble  behavior  and  sound  scatter  at  the  sea  surface,  extended  by  laboratory 
eaqperiments  and  in-situ  measurements  of  near-surface  bubbles  have,  at  last, 
made  it  possible  to  give  reasonably  accurate  predictions  of  surface  duct 
propagation  from  a knowledge  of  wind  speed  over  the  surface  and  duct  velocity 
profile. 

We  will  start  with  a review  of  the  current  oceanographic  description 
of  the  sea  (Section  II),  and  surface  scattering  theory  (Section  III),  supple- 
menting the  basic  formalism  by  the  results  of  recent  laboratory  studies.  The 
theory  of  sound  dispersion,  scatter  and  absorption  by  bubbles,  simple  bubble 
kinetics,  the  extrapolated  results  of  in-situ  near-surface  bubble  measurements 
are  developed  in  Section  IV. 

(U)  Finally  (Section  V),  and  vn.th  the  ver\r  generous  cooperation  of  Dr. 

Homer  R.  Bucker  of  the  Naval  Research  Laboratory,  who  used  his  computerized 
wave  guide  solution^  as  a test  for  our  surface  loss  and  refraction  predictions, 
we  look  again  at  the  data  of  two  SQS-26  propagation  experiments  performed  by 
the  Key  West  Test  and  Evaluation  Detachment.  Bucker 's  wave  theory  solution 
is  modified  by  using  the  observed  wind  speed  at  the  time  of  the  experiment 
to  correct  the  velocity  gradient  for  bubble  presence  and  to  determine  the 
bubble  and  surface  losses  which  now  replace  an  empirical  one  db  surface  loss 
factor  that  he  used  to  obtain  excellent  data  fits. 


II,  THiS  SEA. 


We  assvune  that  the  sea  surface  from  which  the  sound  is  scattered  is 
homogeneous,  and  that  the  temporal  variation  of  sea  heights,  ^ , at  a point 
is  Gaussian  distributed  with  mean  value  zero. 


<c>  = 0 


and  with  variance , a , 


<C'>  - 


The  surface  height  correlation  function  is 


C = “5  = 


(0.0,0)  C(u.v.T) 


where  u,  v,  'T  are  the  spatial  and  temporal  lags. 


The  Fourier  transform  of  2 yields  the  tiiree -dimensional  spectrum  of  the  sea 

^ \ 2 (u,v,T)  exp  [i  (x  u + X V - n t)]  dudvdr  (2) 

*T'^u,v  ^ y 


where  ’^x'^y  ~ components  of  surface  wave  propagation  constant 
O = angular  frequency  of  surface  wave 

we  shall  assume  the  gravity  wave  relation  between  angular  frequency,  propagation 
constant  and  acceleration  of  gravity,  g: 


n = xg, 


J 


This  permits  the  directional  spectrum  of  the  sea  to  be  specified  either 
in  the  two-dimesnsional  rectangular  coordinate  form 


'I'  (x  ,a)  6 (v^-  n) 


(3a) 


or  the  polar  fom. 


X = 6(Vgx-n) 


(3b) 


where  the  Dirac  delta  function,  ^ (Vgx  - Q)  , simply  selects  the  wave  that 
fulfills  the  gravity  wave  relation  between  frequency  and  propagation  constant. 
The  more  commonly  measured  frequency  spectrur.i  is  now  obtained  by 

using  the  energy  density  coordinate  transforaaticn 

2ir 

/ •^j^dKda  = $(J2)  dQ 

0 

This  permits  the  frequency  spectrum  to  be  expressed  as  the  integral  of  the 
polar  form  of  the  X spectrun  over  all  asLiiuth  angles, 

.3  .27T 


$(f2)  r a) 

9 


da 


. (U) 


The  semi-empirical  form  of  the  frequency  spectrum  for  the  fully- 

2 

developed  sea  that  is  most  widely  accepted  today  (Pierson-Moskowitz) 
is 


4»(0)  = 


2 

CC2_ 


exp 


(5) 


I 


where 


= 8.1 X 10 


= Wind  speed  at  19. 5M  above  sea  surface 
= acceleration  of  gravity  in  con-sistent  units 


The  three  statistical  parameters  of  the  sea  surface  that  are  most  significant 

in  specifying  the  sound  scatter  are  the  mean  square  suface  height,  o-^,  the  mean 

2 

square  surface  slope, E , and  the  surface  heirht  correlation  function,  C, 

For  a locally-generated  sea  we  obtain  the  mean  square  height  by 

Xii  voglti  vllig  A r-cvi'-ttii'-jr  op-cv  Vi  Utt  xvi  viic  i aj.  J-jr  — v c oce.^  ^ 


2 ” 

<3^  = ^ 4>(0) 


For  the  mean  square  slope,  it  is  convenient  to  use  the  relation 

3 

determined  from  optical  measurements  of  the  sea 


= 5. 12  X lO’^W.,  + 0.003  + 0.004 
41  ” 


where  W^^is  measured  in  cm/sec  at  a height  of  Ul  feet  over  the  surface. 
Unfortunately  the  measured  speed  is  sensitiv-e  to  the  height  of  the  anemometer; 

u 

however  corrections  can  oe  made. 
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Finally,  for  simplicity  and  because  of  the  gener^  lack  of  two-  | 

dimensional  sea  correlation  data,  we  will  assume  that  the  spatial  correlation  function  | 
is  of  Gaussian  form  and  isotropic: 

C(D  = (8) 

where  I 

6 T ^ I 

X = >U  + V I 

and  L = surface  correlation  length. 

The  assumption  of  an  isotropic  Gaussian  function  makes  it  possible  to  calculate 
the  correlation  length  from  a knowledge  of  rms  slope,  £ , and  rms  height,  o , 
by  the  relation  that  is  valid  for  seas  of  lesser  roughness  ; 


E = V2  ^ (9) 


T 


1 


m.  SOUKD  3GATTiuRIN:i  rHSOKY  LABORATORY  STUDIES 


The  theoretical  solution  to  the  sea  surface  sound  scattering  problem 
usually  starts  with  the  Helmholtz  Integral: 


" 4ir  is  L ‘‘2  an  ■ P Sn  V ^2  /. 


dS 


(10) 


where  p, 


k 
X 

^2 
n 

See  Figure  1, 


the  scattered  acoustic  pressure  at  an  interior  field  point,  Q, 
of  an  enclosed  volume. 

integrating  surface  selected  so  that  acoustic  field  is  zero 
everywhere  except  on  the  erjsonified  water  surface. 

2w/X  = acoustic  propagation  constant 
acoustic  wave  length  in  water 

distance  from  surface  scattering  region  to  field  porition. 


= no 


rrriQl  to  surfoc^  ^ 


e scattering  point. 


After  applying  the  boundary  conditions  for  acoustic  pressure  and 
acoustic  particle  velocity  at  the  sea  surface,  and  assuming  that  there  is  no 
shadowing  or  secondary  scattering  even  for  near-grazing  incidence,  we 
obtain  the  integral  for  the  scattered  pressure  (now  called  AP2  r Pq)  from  a given 

surface  of  scattering  area  A A in  terms  of  the  geometry: 


plane 


(11) 


where  A p. 


'plane 


scattered  pressure  for  a pl3j;s  surface  A A ov^r  which 
there  is  negligible  curvature  of  the  incident  wave  front, 


= ikp,  AA  cos  9,  e^^^'^2  ^i^Vz^rr, 


AA 


= insoiilfied  STirface  area  over  which  is  constant  and  tlJB 
sagitta  of  the  wave  front  is  less  than  A/8  . 


F 


= incident  sound  pressure  at  sea  surface. 

= incident  sound  angular  frequency 

1 + cos  0^  cos  02  ” sin  0^  sin  02  cos  0 
cos  0j  (cos  0j  + cos  02) 


0 

6 

0 


1 

2 

3 


= angle  of  incidence  (measured  with  normal) 

= angle  of  scatter  (measured  with  normal) 

= azimuthal  angle  of  scatter  plane  with  respect  to 
incident  plane. 

= k (ij  sin  0^  - cos  0^) 

= k (ij^  sin  02  cos  63  + 12 


•^1  '^2  '-3 


are  the  unit  normals  in  the  x,  y and  z (depth)  directions, 
respectively. 


It  should  be  noted  that  the  solution  being  presented  assur.es  that  the  illu;ninatcd 
surface  is  in  a snail  section  {/Ka  <<  rj_)  of  the  Fraunhofer  legion  of  the  source 
field  and  that  the  field  position  is  in  the  Fraunhofer  region  of  the  illuninated 
area  (1*2 » r2  >>  AA/x  ),  This  point  is  emphasized  by  our  A notation  for  the 
illVLidnated  area. 

The  mean  square  scattered  pressure,  which  in  general  includes  both 
coherent  and  incoherent  sound  intensity,  is  formulated: 


(AA)  . . 

X y x'y 


(dxdx'dydy') 


8 


r 


where 


K 


-1  ■ ~2 


u 


I 


pi 

* 

< > 


X - X 

y'  - y 

Ap 


Ap  * 

2 ^2 
plane  plane 


complex  conjugate 
ensemble  average 


The  ensemble  average,  < ) , can  be  easily  evaluated  for  our  surface 
with  its  two-dimensional  Gaussian  distribution: 


^ giK^,  (C  - C ') ^ ^ [kCT  (cos  6^  + cos  6^)  ^ Cl-C] 


Define  the  ACOUSTICAL  ROUGHNESS  OF  THE  SEA  SURFACE  = g 


g = [ko  (cos  + cos  6^)  ] “ 


(12) 


so  that  the  mean  scattered  intensity  can  be  written 


/A  A T F ^-g(l'C) 

<Ap  Ap  *)  = I - e 

^ ^ (AA) 


$ S S * V*  dxdx'dydy' 


X y X y 


IL 


If  our  interests  were  to  be  restricted  to  low  roughness  (g  < 1)  at  this 
point,  it  would  be  useful  to  use  the  approxLnation  e — (e“S)(l+gC), 

This  would  pemit  inraediate  separation  into  a coherent  tern,  e”S,  and  an 
incoherent  term  which  the  integration  would  transform  into  the  two  dimensional 
surface  spectral  description,  eqn.  ("b)  » (Aciner-Khinchinc  ti'ansformation) » 
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The  rarity  of  experimental  measurements  of  two-dimensional  surface  spectral 

descriptions  turns  us  away  from  this  approach. 

Instead,  for  scatter  in  the  specular  direction,  it  is  convenient  to 

try  a series  approximation  to  the  assumed  isotrooic  Gaussian  spatial  correlcition 

5,6 

function  (Eq.  8).  The  inte<p7ation  then  produces  the  relative  scattered 
intensi^’y,  (compared  with  the  plane  surface  reflected  intensity)  specialized 
to  speculax'  scatter,  in  which  6^  = 6^  . 0^  = 0 » F = 1 , g = iik^2 
where  = s grazing  angle. 


<pp*> 


(Ap^Ap^*) 

^plane 


(13) 


The  two  terras  of  Eq. (13)  represent  coherent  scattered  sound  (the  e~S 
terra)  and  incoherent  scattering  (the  second  terra).  ^ study  of  the  terra 
in  parentheses,  defined  as 


S(g) 


ra 


(li^) 


shows  that  S (g)  •>  g for  g « ! 

and  S (g ) ■>  ~ for  g » 1 

2 

Since  it  can  be  assumed  that  L < AA,  we  can  conclude  that,  in  the  limits, 
the  specular  scatter  for  a Gaussian  sea  will  be: 


For  Small  Acoustical  Roughness,  g <<  1:  (pp*>  = e 

nL^  IT 

For  Large  Acoustical  Roughness , g » 1:  {pp*)  = -r-r  = — t—z 

E AA  , » . 2^ 

„k  E AA.  sin  (f: 

In  practice  the  coherent  tern  is  dominant  ’pproxLraately  up  to  g-1,  anc 

sometimes  beyond  that  roughness,  depending  on  the  beam  widt'  . Fig.  2 

10 


(15) 

(16) 


shows  the  relative  behavior  as  a function  of  acoustical  ro'.ighness. 

It  is  seen  that  for  snail  acoustical  roughness,  the  intensity  of  specular 
. scatter  is  the  factor  e“^  tines  the  mirror  reflection  , and  the  concept 

I 

of  a reflcctio:j  coefficient  is  useful. 

However,  it  is  important  to  observe  also  that  the  value  of  dB  loss  per 

bounce  is  in  general  an  inadequate,  and  in  fact  ;.;isleading,  concept  for  the 

rough  surface  case,  when  g>^  1 the  specific  value  of  the  scattering  area 

enters  into  the  as^Taptotic  value  of  (pp*)  wnich  will  therefore  chaiige  from 

one  bounce  to  the  next  as  the  sound  wave  continues  to  diverge  and  scatter, 

when  the  iYaunltofer  conditions  are  violated  (large  beaiii  width)  it  is 

necessary  either  a)  to  aad  the  scatter  contributions  from  the  several  areas 

A A or  b)  to  re-sol ve  the  problem,  assaying  a diverging  ;.ave  at  the 

T 0 

very  beginning  of  the  derivation.  * Ihe  latter  attack  and  its  verification 

Q 

by  laboratory  model^  show  that  for  specular  scatter  from  an 

exponentially  correlated  surface  there  is  virtually'-  no  difference  between  the 

Fraunhofer  and  iYesnel  solutions  if  g<l,  i;owe\'cr,  for  g>l  the  scattered 

inten.sity  may  'oe  significantly  greater  than  predicted  by  the  single  Fraunliofer 

solution,  T.;is  is  cecause,  com.ionly,  significant  signal  v;ill  e:d.st  ove.r  an  area 

containing  more  than  one  A A as  defined  in  our  deriwition. 

; Laboratory  models  have  been  used  in  order  to  supplement  the  capabilities 

9 10 

of  current  theory,  *'  We  have  recently  studied  the  behavior  of  the  coherent 
component  for  cross-wind,  up  and  down  wind  specular  scatter  under  widely-varying 
laboratory  seas.  Our  "sea.®"  are  produced  by  combinations  of  a paddle  and  five 
fans  blowing  air  over  the  surface  of  our  25'  x 6'  x 8'  deep  anechoic  tank. 

In  this  way  we  obtain  water  waves  that  give  us  scaled  replicas  of  canpound  sea 
and  swell  surfaces  that  are  realistically  time-varying,  anisotropic,  and  near- 
Gaussian.  The  probability  density  function,  spatial  and  tern  >oral  correlations 

' ■"The  Gaussian  correlated  surface  would  give  similar  res'olts. 


11 


and  frequency  spectrum  ccanpare  closely  to  ocean  measurements  , except  for 

scale  factors  for  frequency,  time  and  space  lag,  of  course. 

Figure  3 shows  three  of  the  laboratory  sea  and  swell  spectra.  The 

probability  density  functions  of  the  wave  height  for  the  three  cases  are  shown 

in  Fig.  li;  the  near-Gaussian  probability  density  function  is  very  close  to 

that  observed  at  sea.^^  Equation  (13)  states  that  the  ratio  of  the  intensity 

of  the  coherent  component  to  the  intensity  reflected  from  a mirror-like 

surface  should  be  e”®.  Figure  5 shows  this  ratio  divided  by  e"®  over  four 

decades  of  acoustical  roughness;  for  g£  1 the  ordinate  is  essentially  unity. 

This  excellent  verification  of  theory,  when  tested  under  widely  varying 

laboratory  seas,  has  encouraged  us  to  apply  it  to  the  real  ocean.  Other 

11  12 

laboratory  model  results  * suggest  that  for  g<l,  at  near  grazing  incidence, 
a shadowing  correction  is  not  necessary. 


IV.  EFFECTS  OF  NEAR-SURFACE  BUBBLES 


Bubbles  near  the  sea  surface  can  modify  the  forward  scatter  of  sound 
In  two  ways:  (a)  By  near-surface  refraction;  the  presence  of  all  bubbles 

(resonant  or  not)  will  create  a below-surface  medium  through  which  the  speed 
of  sound  differs  from  that  in  the  non-bubbly  water;  this  change  of  speed  will 
refract  incoming  rays,  modify  the  angle  of  incidence  at  the  rough-scattering 
surface  and  change  the  distance  betweeen  bounces  in  channel  propagation. 

(b)  By  resonant  absorption  and  scatter;  resonant  near-surface  bubbles 
can  cause  attenuation  that  is  generally  unresolvable  from  rough  surface 
scatter  loss. 

The  objectives  of  the  present  work  are  arbitrarily  restricted  to 
"surface'*  losses  at  frequencies  below  16  kHz.  Nevertheless,  because  bubble 
refraction  depends  on  the  total  volume  of  air  in  bubble  form,  it  is  necessary 
to  consider  the  numbers  of  near-surface  bubbles  of  all  radii. 

In  soite  of  the  creat  interest  in  the  subject  for  almost  30  years, 
the  literature  of  ocean  bubble  data  is  very  thin:  In-situ  measurements  of 

bubbles  in  off-shore  brea'<ing  waves  have  been  reported^ bubble  generation 
and  persistance  by  an  artifical  wind  over  a laboratory  tank  filled  with  sea 
water  has  been  described^^*;  inferred  bubble  densities  have  bet-  calculated 

l5 

from  certain  AKCS  data  ; in-situ  acoustical  measurements  in  60  ft  of  iso- 
thermal water  in  Mission  Bay,  San  Diego  have  been  summarized^  j a laboratory 

study  of  surface  bubbles  has  yielded  data  produced  by  colliding  waves  of 
17 

sea  watero  It  will  be  necessary  for  us  to  extrapolate  beyond  these  bare 
measurements  available. 

A,  Dependence  on  Radius 

Because  of  experimental  difficulties,  in  no  case  are  there  published 
in-situ  measurements  of  ocean  bubbles  of  radius  greater  thin  200 

13 
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microns  (resonant  frequency  less  than  approximately  15  kHz)  or  less  than 
roughly  20  microns  (resonant  frequency  greater  than  150  kHz) . Most  of  the 
available  data  fit  approximately  the  solid  line  in  Fig.  6 > plotted  on 
the  basis  of  fractional  volume  of  air  in  bubble  form  to  water  u(R)dR.  From 
the  fractional  volume  it  is  easy  to  calculate  the  number  of  bubbles  per  unit 
volume,  in  radius  increment  dR,  by  using 


n(R)  dR 


u(R)  dR 
4/3ttR^ 


(17) 


There  is  some  evidence  that  the  bubble  populations  break  up  roughly  into  two 
groups:  those  smaller  than  approximately  60  micron  radius  postulated  as 

attached  to  motes  and  due  principally  to  biological  activity,  cosmic  rays, 
captured  aerosols  and  photosynthesis  of  phytoplankton,  and  the  larger  bubbles 
which  we  assume  come  from  breaking  waves  (and  decaying  matter  on  the  sea 
floor  in  shallow  water) . 

V.  To  calculate  the  near-surface  refraction  due  to  bubbles  we  have  had  to 
extrapolate  (dashed  lines)  the  data-based  curves  of  Fig,  6,  which  were 
measured  with  incipient  breaking  waves  at  wind  speeds  U-6  knots  at  3*3  meter 

depth.  For  radii  smaller  than  30a  we  extrapolate  to  one  micron  radius  at  the 
-1/2 

K2R  ' dependence  that  has  been  found  for  the  smallest  measured  ocean 
bubbles^^  and  in  tap  water^^j  below  1 micron,  at  slope  following  behavior 

19  20 

loosely  suggested  by  laboratory  cavitation  studies,  * For  bubbles  larger  than 
60u  we  have  used  the  dependence  observed^  to  100  microns  and  have 


Hi 


-3 

extrapolated  to  largest  radii  at  Blanchard  and  Woodcock's  slope  of  K^R 
The  radius  at  which  one  shifts  to  the  suggested  slope  for  largest  bubbles 
is,  admittedly,  somewhat  arbitrary.  That  slope  starts  at  about  50  microns 
in  the  breaking  wave  data;  however,  a re-analysis  of  Glotov's  data  shows  that 

_3 

the  K^R  behavior  is  confirmed  for  R > 120u  . The  constants  for  the  slopes 
of  Fig,  6 are  shown  in  Table  1. 

TABLE  ONE 

Estimated  Fractional  Bubble  Air  to  Water  at  3.3M  depth  for  Winds  Approximately 


\ 

6 Knots 

Radius  Range 
Microns 

Dependence 

u(Rq.3.3M) 

Constant 

[dimensions] 

/ 

dR 

^0 

1.0 

- k 0 

'^I'^O 

Ki  = 

1.5  X 

Iq-9 

1 J 

1.5  V 

10"' 

1. 

0 <, 

^0 

<_  30. 

K R 

*^2  = 

1.5  X 

10-' 

13.8  X 

10"' 

30. 

< 

R 

0 

£ 60. 

8.5  X 

10"' 

[y°3 

5.8  X 

10"' 

60. 

< 

*^0 

<_  100. 

Vo*' 

^4  = 

2.5  X 

lO"^^ 

[u"^] 

7.7  X 

10"' 

100. 

< 

^0 

Vo"' 

S = 

2.5  X 

10"^ 

12.4  X 

10"' 

= A.l  X 10 


Total  Effective  U 


B.  Dependence  on  Depth 


For  our  purposes,  interest  lies  in  the  bubble  densities  (and  fractions) 

as  a function  of  depth  at  a fixed  frequency,  rather  than  for  fixed  radius. 

Studies  in  well-mixed  ocean  water  suggest  that  for  bubbles  larger 

than  approximately  30  microns  (sea  level  resonant  frequency  less  than  100  kHz) 

buoyant  forces  and  turbulent  entrainment  determine  the  bubble  life  time 

rather  than  gas  diffusion  through  the  bubble  walls.  If  we  formally  make  the 

simplifying  assumption  that  these  bubbles  can  be  lost  only  at  the  ocean  surface, 

then  bubbles  of  surface  (3.3  Meter  depth)  resonant  frequency  f^,  radius 

will  be  found  at  depth  Z in  the  same  numbers,  but  at  the  isothermally*^ 

1/3 

compressed  radius,  = Rq6  and  resonant  frequency  f^^  where 

P “ 1.33/(1  + O.IZ). 

For  resonant  air  bubbles  kR^  = ~ where  c « 1 = 5 x lo^  c.-Aec 

y = 1.4,  ? = 10^  ^1+Z/lO  dyne/cm'^,  °p^  » 1.0  gm/cm^,  Z = depth  in 

meters , 

The  bubble  fraction  at  depth  Z will  be 

u(Rqz,Z)  = 


Analytically,  our  assumption  for  large  bubbles  is:  from  the  surface 
to  depths  where  the  entrained  bubble  drag  force  is  equal  to  or  greater  than 
the  buoyant  force 

n(R^,3.3M)  dR  = n(R^2»2)  dR ^ 
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80  that  following  isothermal  compression, 


u(Rq2,Z)  = [u(Rq,3.3M)][B]. 


These  bubbles,  however,  have  a resonant  frequency  which  is  not  but, 
because  of  changed  radius  and  pressure,  is 


^OZ  “ ^0  ^ 


We  can  now  calculate  the  gas  fraction  of  bubbles  at  the  desired 
resonant  frequency  at  depth  Z,  u(f^,Z)  , by  assuming  that  the  frequency 
dependence  within  a given  band  of  bubble  sizes  is  the  same  as  at  the  surface, 
chat  is, 


u(f  ,3.3H) 

0 


k „ f ‘ 
no  u 


u(f,-.^rZ) 


'‘n2  'oz 


where  n is  the  frequency  dependence  within  a particular  frequency  band. 


Therefore 


u(fQ.Z)  - 


Finally,  from  ( l8)  and  (19) 


u(£„.Z)  - u(f„.3.3M) 


and,  taking  account  of  the  change  in  radius, 

„(f^,Z)  - n(f(,.3-3M) 


liLi. 


The  predicted  depth  dependences  for  the  three  most  important  bands 


of  surface  radii  are  shown  in  Table  2, 

TABLE  2 

Fractional  Air  in  Bubble  Form  and  Number  of  Bubbles  per  Unit  Volume 
as  a Function  of  Depth  for  Large  Bubbles 


^0 

^0 

n 

u(fQ,Z)/u(fQ,3.3M) 

N(fQ,Z)/N(fQ,3.3M) 

microns 

kHz 

_ g5/6  n+1 

-5/6  n + 5/2 

■ P 

30  to  60 

60  to  100 

+1 

gll/6 

glO/3 

60  to  100 

32  to  60 

-1 

3^/3 

> 100 

less  than  32  +3 

3^/2 

3^ 

It  is 

clear  that  the 

assumptions 

that  lead  to  the  values 

in  Table  2 

can  only  be 

considered  as  a 

perturbation 

analysis  of  the  change 

of  bubble 

populations 

with  depth  since 

they  imply 

inconsistencies,  particularly  at  the 

edges  of  the 

bands,  as  depth 

is  increased.  Nevertheless,  the  only  bubble 

density  data 

available 

plotted  against  the  absolute  pressure,  1.333 

for  four  frequencies  between  40  and  96  kHz  show  slopes  of  approximately  6 , 

1.7  3.3 

in  reasonable  agreement  with  the  predictions  of  6 " to  6 ’ . 

When  we  plot  the  frequency  band  contributions  to  the  total  bubble  fraction, 
( / u R dr),  as  a function  of  total  pressure  parameter  6 the  total  bubble 

fraction  U(z)  at  wind  speed  6 knots  is  found  to  vary  (with  less  than  deviation) 

+2  +1 
as  6 from  the  surface  to  approximately  lOM  depth,  blending  into  3 

at  greater  depths  (near  20  Meters) . Therefore  the  dpeth  dependence  that  we 

will  use  for  the  near-surface  refraction  calculation  is: 

U(Z)  = 4.1  X lO"®  3'*’^  t22) 

. ■ «*  ^ 


that  the  number  of  bubbles  shows  "an  abrupt 


) 

I 


C.  Dependence  on  Wind  Speed 

It  has  been  found  ^ 

Increase  ...  approximately  ...  exponential"  for  winds  above  about  6 knots. 

21 

The  comment  agrees  roughly  with  the  curve  of  foam  coverage  vs.  wind  speed 
which  can  be  written  approximately  as: 

U(Z=0,W)  = U(Z=0,6  knots)  [1  + 10"^(W/10)^  (23) 

where  W = wind  speed  in  knots. 

There  is  a second  wind  effect:  the  larger,  wind-generated  bubbles  will  be 

brought  to  depths  by  increased  turbulent  entrainments  A complete  knowledge  of 

this  behavior  may  be  possible  when  there  is  a better  understanding  of  under- 

22 

water  motion  as  a function  of  wind  speed,  insolation  and  depth. 

However,  the  factor  of  bubble  entrainment  by  water  motion  can  be  crudely 
evaluated  by  the  ioilowing  simple  argument: 

A bubble  will  be  drawn  down  from  the  surface  when  the  drag  force  is 
greater  than  the  buoyant  force. 

The  convective  velocity  in  the  medium  will  probably  be  Gaussian  distributed. 
Nevertheless,  the  mean  downward  con'/ective  velocity  will  increase  with  wind 

23  2J, 

speed  as  does  the  velocity  spectrum  and  the  depth  of  the  well-mixed  layer. 

We  assume  that  the  mean  convective  velocity  is  prcporticnal  to  the  wind 
speed  and  that  the  mean  number  of  bubbles  drawn  downward  under  a given 
condition  of  underwater  turbulence  is  proportional  to  the  mean  velocity;  the 
dependence  on  radius  will  be  ignored  at  this  time.  This  "zero  order" 
estimation  of  the  effect  of  turbulent  entrainment,  combined  with  eq.(23) 
for  the  increased  foam  available  leads  to  the  wind  dependence: 


U(Z,Rq,W) 
U(Z,Rq,6  knots) 


»(W/6)  (1  + io'iw/lo)^j 


(2k) 
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D.  Summary  of  Bubble  Dependences 

Combining  the  dependences  on  radius  of  bubble,  bubble  depth  and  wind 
speed,  we  get  the  total  bubble  fraction: 

U(Z,W)  » [i4.1  X lO"®^^]  [W/6J  [1  + 10'^(W/10)^]  (25) 

It  is  appropriate  to  point  out  that  other  bubble  factors  exist  in 
addition  to  the  effects  that  we  have  considered.  Two  significant  elements 
are  the  dependence  on  dayli^t  and  precipitation;  others  are  dependence  on 
water  temperature  and  gradient,  percentage  gas  saturation,  marine  biological 
activity.  As  examples:  Blanchard  and  Woodcock  have  shown  that  snow  and 

rain  produce  very  large  numbers  of  bubbles.  Their  penetration  in  laboratory 
experiments  was  no  more  than  2 cm  below  the  water  surface  but  it  would  be 
expected  that  these  bubbles  would  penetrate  much  farther  in  the  presence  of 
convection  cells  accompanying  storms  at  sea.  Small  bubbles  (less  than  about 
60V  ) are  about  twice  as  populous  in  the  daytiii.e  as  at  night.  It  has  been 
postulated  that  the  increased  numbers  may  be  due  to  photosynthetic  action. 
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E.  Theory  of  Dispersion  and  Absorption  by  Bubbles 

Following  classical  developments  tin  incident  acoustic  pressure, 
of  wave  length  very  large  compared  to  the  bubble  radius  R,  having  a velocity 
potential  will  produce  a total  field 

N 

(ir)  = Ip.  (ir)  + X ("^)  "^n*  (26) 

^ n = 1 ' ‘ ' 


where  ir  = coordinate  position  of  scattering  region 

ir-ir  = distance  between  bubble  and  center  of  scattering  region 
n 

k =»  2n/\  = propagation  constant  in  bubble-free  water 

Rip 

A = — 

^ {‘*’0/")^  - 1 - i6 

R = bubble  radius 

w = incident  sound  frequency 

= resonant  frequency  of  bubble 

6 = damping  constant  of  bubble 

The  average  over  all  possible  sizes  and  configurations  of  bubbles  is 

"tr'-  ir  'r 

where 


Equation  (2?)is  the  integral  equation  formulation  of  the  differential 
equation 


{il>)  + ill))  = -4ffG<0>  or 

<0>  + (r)  <0>  = 0 


where 


k|,(r) 


n fr.R)  R d R 
(a,o/o,)2.  i-i6(R) 


(28) 


The  quantity,  — r — = 1 + — ^ , is  the  index  of  refraction  in  the  bubbly 
region.  Because  of  the  presence  of  the  damping  term,  the  index  will  have 
an  imaginary  part  which  comprises  the  attenuation  due  to  absorption  and 
scatter  in  the  bubbly  water  in  addition  to  the  real  part,  defining  the  local 
speed  of  sound.  These  terms  are: 


Re:  |k,  (r)|  = k + 


2tt 


/ 


R n(R) 


ll  d R 


[(^)^  - 


+ 6' 


(2?) 


Im:  lyr)}  = f J 


6 R n (R)  d R 


(30) 


F.  Predictions  of  Refraction  by  Bubbles 

The  dispersion  relation  Eq.  (29) » since  it  predicts  a decrease 

of  sound  speed  caused  by  bubbles  of  resonant  frequency  greater  than  the 

Incident  frequency  and  an  increase  of  speed  for  bubbles  with  lower  resonant 

frequencies,  would  require  precise  information  of  bubble  populations  for 

an  evaluation  for  all  sound  frequencies.  However,  if  our  interest  is  restricted 

to  sound  frequencies  less  than  16  kHz,  integration  of  Fig.  6 shows  that 

bubbles  under  that  resonant  radius  (R^  < 200.  microns)  represent 

approximately  93%  of  the  total  fractional  air  volume.  For  the  calculation 

of  refraction  at  frequencies  less  than  16  kHz  we  consider  only  the  effect  of 

bubbles  of  resonant  frequency  freater  than  the  driving  frequency. 

2 

When  >>  w , we  drop  6 in  the  denominator  of  Eq.  (29l,use 
definition  of  u (R)  to  obtain  the  simplified  form: 

Re: 

Since  the  expression  in  parentheses  is  no-rly  unity,  v?e  find 
c^  = c 1 - 8.0  X 10^  U (1  + Z/10)"^ 


(31) 


(32) 


where  U(Z,U)  is  given  by  equation  (2^). 
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The  significance  of  this  departure  of  the  speed  from  the  bubble-free 


value  is  evident  when  we  insert  Equation  (25)  and  look  at  the  gradient 

of  the  near-surface  propagation  speed.  For  the  reference  wind  speed  of 

-4 

6 knots  the  gradient  of  the  velocity,  proportioanl  to  (1  + O.IZ)  , has  a 
surface  value  of  0.26  sec  decreasing  to  0.016  sec  ^ by  10  Meter  depth. 

Since  the  velocity  gradient  due  to  pressure  in  isothermal  water  is  0.018  sec  ^ 
our  conclusion  must  be  that  ray  curvature  in  approximately  the  top  10  meters 
is  very  strongly  influenced  by  bubbles  even  for  winds  as  low  as  6 knots. 

Because  of  bubbles,  rays  approaching  the  surface  will  refract  upward, 
meeting  the  surface  in  a cusplike  ray  path,  and  making  more  bounces  for  a 
given  horizontal  range  than  for  a bubble-free  region. 


which  can  he  re-expressed  as 

3 - 1.50  X 10^  (1+2/10)“^'^^  u(fQ,Z,W)  db/Meter  (3U) 

b 

The  conclusion  of  the  calculation  of  "surface  loss  in  db/bounce 
due  to  bubbles  requires  integration  of  Eq.  (3U  ) over  the  surface  duct  ray 
path  of  varying  curvature  and  varying  attenuation: 


V.  SURFACE  DUCT  PROPAGATION  - COM  ARISON  WITH  EXPEItlMEin’ 


1 

Sucker  has  recently  solved  the  surface  duct  propagation  problem  by 
a wave  theory  approach.  His  modal  components,  which  correspond  to  up  and 
down  propagating  waves,  with  appropriate  ray  directions,  lose  energy  when 
they  interact  v?th  the  surface  and  by  leakage  through  the  thermocline.  Our 
proposal  here  is  to  introduce  the  results  of  our  development  in  Sections  III 
and  IV  int<^  Sucker's  sclut7  -n.  Two  changes  are  necessary’':  the  velocity 

gradient  in  the  duct  will  now  be  due  to  bubbles  as  well  as  temperature  and 
pressure  effects;  the  loss  at  the  surface  will  l>e  made  up  of  the  losses  per 
bounce  due  to  below-surface  extinction  of  sound  by  bubbles  and  simple,  at-surface, 
coherent  specular  scatter  (because,  in  gen'ral,  g<l),  calculated  for  each 
contributing  mode.  The  surface  "reflection  coefficient"  is  thereby  directly 
attributable  to  the  wind  velocity  which  is  the  sole  source  of  the  rms  v;ave 
height'*"  in  the  examples  chosen  and  a principal  cause  of  the  bubbles. 

First,  let  us  lock  at  how  the  losses  vary  with  wind  soeed  and  ray  angle 

of  incidence.  Our  loss  calculations  are  based  on  eqs.  l5,  ;5>  36.  Our  ray 

paths  are  calculated  from  the  bubble-caused  £yadient  using  cqs.  25,  32  together 

with  a temperature  plus  pressure-caused  gradient  of  OoO!48  i Fig-are  7 

illustrates  the  relative  contributions  of  surface  loss  and  integrated  bubble 

loss  for  a ranee  of  winds  with  a ray  that  is  horizontal  at  UO  M,  *?  = U.38°. 
belCT.v 

The  frequency  is'  16  kHz.  It  is  observed  that  the  bubble  contribution 
dominates  only  for  low  wind  speeds. 

*In  the  data  report^'^  the  "seaman's  eye"  estimates  of  the  neak-to-peak  wave 
height,  H,  ar^  in  good  (BQX)  arreenent  with  the  calculatic:-.  from  equation  (6) 
when  it  is  assumed  that  H ■*  60-  , for  Gaussian  wave  systems.  However,  our 
laboratory  experiments  suggest  that,  in  general,  the  total  height  due  to 
swell  and  "sea"  should  be  used  to  determine  tr. 


1 


Fig^e  8 shows  that  the  bubble  loss  per  bounce  is  not  very  sensitive 
to  the  depth  of  the  horizontal  ray  (or  angle  of  incidence  at  the  surface); 
the  loss  per  bounce  at  the  surface,  however,  increases  as  the  square  of  the 
sine  of  the  grazing  angle. 

For  higher  frequencies  it  would  be  expected  that  the  neglect  of  the 
incoherent  surface  scattering  term  would  no  longer  be  possible;  the  forvrard 
scatter  loss  would  reach  an  asymptotic  value  which  vrould  be  smaller  than  the  increased 
bubble  loss.  For  frequencies  lower  than  our  exajnple  the  coherent  scatter 
approximation,  e”^,  becomes  even  better,  and  would  be  applicable  to  hi^er 
sea  states;  the  bubble  loss  would  be  expected  to  continue  to  be  smaller  than 
the  surface  loss. 

The  exoerlmental  propagation  data  for  our  study  are  again  from  Runs  73 
and  99  of  the  oper~tional  tests  of  the  SQ.S-26  sonar^"^  in  which  both  source 
and  receiver  wex'e  in  the  ciuet.  The  wave  heights  a:  e such  tluio  g < i lor  an 
effective  modes  of  the  duct  solution. 

Since  Sucker's  duct  solution  assumes  constant  velocity  gradients,  we 
have  had  to  approximate  the  calcuLated  total  gradient,  which  goes  as  ^ , 

by  splitting  it  into  constant  sections  and  G2  iFig,  9).  The  near-surface 
section,  G]_  , has  been  given  a value  twice  the  uncorrected  gradient  due  to 
temperature  and  pressure  and  the  remainder,  G2  , has  been  assumed  to  have  no 
bubble  influence.  The  values  are  given  in  Table  3»  In  tV;is  '..'ay,  the  duct 
propagatiun  loss  curves  have  been  recalculated  bv  Bucker.  Thev  compare  well 
with  the  data  (Figs.  10,  11 ).  Ihu  uashod  lines  in  Jig-ures  10  and  11  are  the 

rruch  ~'r)orcr  prodicticns  based  on  the  equations  of  "ro-ject  AX~  / 
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TABLE  3 


Data  for  Two  Surface  Duct  Experiments 


Run  No. 

Vs 

Zt 

^2 

W 

M/sec 

M 

sec“^ 

sec~l 

knots 

73 

1513.1 

36.6 

0.096 

0.0U8 

17 

99 

1510.9 

51.8 

0.102 

0.03U 

15 

where 

^s 

= velocity 

at  surface 

Zt 

■ depth  of 

thermocline 

= velocity 

gradient,  surface  to  8 yards 

^2 

= velocity 

gradient,  8 

yards  to 

thermocline 

Other,  higher  values  of  the  assumed  bilinear  gradient  and  a closer- 
fitting tri-linear  gradient  have  also  been  consideired.  However,  the  resulting 
propagation  prediction  curves  show  excessive  losses  compared  to  exjjerimento 
Since,  at  these  wind  speeds,  the  bubble  less  is  less  than  the  surface  loss, 
the  velocity  gradient  can  be  a major  influence.  This  suggests  that  O’ur  total 
bubble  populations,  U(fQ,Z,V/),  may  be  too  high,  perhaps  by  a factor  of  two, 
possibly  due  to  the  wrong  wind  speed  dependence  having  been  assumed. 


VI . CONCLUSION 


The  purpose  of  this  paper  has  been  to  show  that  for  a large  range  of 
ocean  surface  roughness,  g < 1,  it  is  possible  to  predict  surface  duct 
propagation  by  using  the  simplest  surface  scattering  theory  (coherent  terra 
only)  and  estimates  of  bubble  presence  as  inputs  to  the  Sucker  modal  solution. 
Nevertheless,  it  is  the  necessity  of  knowing  more  about  the  ocean 
parameters  of  sound  propagation  that  has  been  demonstrated  by  this  paper. 

It  is  hoped  that  the  present  work  will  ser/e  as  a spur  for  more  fundamental 
studies  to  replace  the  speculative  assumptions  that  we  have  been  forced  to 
make  about  bubble  dependence  on  wind  speed. 
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